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Abstract We review here briefly the nature of precursors of global failures in three 
different kinds of many-body dynamical systems. First, we consider the 
lattice models of self-organised criticality in sandpiles and investigate 
numerically the effect of pulsed perturbations to the systems prior to 
reaching their respective critical points. We consider next, the random 
strength fiber bundle models, under global load sharing approximation, 
and derive analytically the partial failure response behavior at loading 
level less than its global failure or critical point. Finally, we consider 
the two-fractal overlap model of earthquake and analyse numerically 
the overlap time series data as one fractal moves over the other with 
uniform velocity. The precursors of global or major failure in all three 
cases are shown to be very well characterized and prominent. 

Keywords: Fracture, earthquake, avalanches in sandpile, self-organised criticallity, 
fiber bundle, fractals. Cantor sets. 



1. Introduction 

A major failure of a solid or a dynamic catastrophe can often be viewed 
as a phase transition from an unbroken or non-chaotic phase to a broken 
or chaotic phase. Some obvious correlations existing in the system before 
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the failure, are lost in the failure process. However, the equivalence can 
be made precise and in fact the critical behavior goes to the bone of any 
failure or catastrophic phenomena. 

Several dynamical models of cooperative failure dynamics [11 are now 
well-studied and their criticality at the global failure point are well es- 
tablished. Here, we have reviewed some of the numerical studies of 
precursor behavior in two self-organising dynamical 1 models as one 
approaches the self-organised critical (SOC) point of the avalanches in 
the lattice sandpile models. The critical behavior of a random fiber 
bundle model of failure or fracture of a solid under global load shar- 
ing is now very precisely demonstrated. The dynamics of failure also 
become critically slow there and all these have been demonstrated ana- 
lytically. The precursors of the global failure in the fiber bundle models 
can therefore be discussed analytically. The resulting universality of the 
surface roughness in any such fracture process has also been well docu- 
mented and analysed. The two-fractal overlap model |H] is a very recent 
modeling approach of earthquake dynamics. The time series of overlap 
magnitudes obtained, when one fractal slides over the other with uni- 
form velocity, represents the model seismic activity variations. This time 
series analysis suggests some precursors of large events in this model. 

2. Precursors in SOC models of sandpile 

A 'pile' of dry sand is a unique example SOC system in nature. A 
growing sandpile gradually comes to a 'quasi-stable' state through its 
self-organising dynamics. This 'quasi-stable' state is called the critical 
state of the system as the system exhibits power law behavior there. 
Because of avalanches of all sizes, this critical point in the pile is also 
a catastrophic one. The dynamics of growing sandpiles are successfully 
modeled by the lattice model |^ of Bak, Tang and Wisenfeld (BTW) 
in 1987. A stochastic version of sand pile model has been introduced 
by Manna [Sj which also shows SOC, although it belongs to a different 
universality class. Both the models have been studied extensively at 
their criticality. Here, we study the sub-critical behavior of both the 
models and look for the precursors of the critical state. 

BTW Model. Let us consider a square lattice of size LxL. At each 
lattice site there is an integer variable hij which represents the 

height of the sand column at that site. A unit of height (one sand grain) 
is added at a randomly chosen site at each time step and the system 
evolves in discrete time. The dynamics starts as soon as any site 
has got a height equal to the threshold value {hth= 4): the site topples, 
i.e., hi^j becomes zero there, and the heights of the four neighboring sites 
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increase by one unit 

hij hij — 4, hi±ij hi±ij + 1, hij±i hi,j±i + 1- (1-1) 

If, due to this toppling at site (i, j), any neighboring site become unstable 
(its height reaches the threshold value), they in turn follow the same 
dynamics. The process continues till all sites become stable (/lij < hm 
for all When toppling occurs at the boundary of the lattice (four 

nearest neighbors are not available), extra heights get off the lattice and 
are removed from the system. 

With a very slow but steady rate of addition of unit height (sand 
grain) at random sites of the lattice, the avalanches get correlated over 
longer and longer ranges and the average height {hav) of the system grows 
with time. Gradually the correlation length (,^) becomes of the order the 
system size L. Here, on average, the additional height units start leaving 
the system as the system approaches toward a critical average height 
hc{L) and the average height hav remains stable there (see Fig. 1(a)). 
Also the system becomes critical here (for L — > oo) as the distributions of 
the avalanche sizes and the corresponding life times follow robust power 
laws 15 . 
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Figure 1. The growth of average height hav against the number of iterations of adding unit 
heights (L = 100). Eventually hav settles at hc{L). In the inset, we show the finite size 
behavior of the critical height hc{L), obtained from simulation results for different L. (a) For 
the BTW model; (b) for the Manna model. 

Here, a finite size scaling fit hc{L) = /ic(oo) + CL^^/'^ (obtained by 
setting ^ ~ I hc{L) — /ic(oo) \^'^= L), where C is a constant, with i/ ~ 
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1.0 gives he = hc{oo) ~ 2.124 (see inset of Fig. 1(a)). Similar finite 
size scaling fit with = 1.0 gave /ic(oo) — 2.124 in earlier large scale 
simulations [7]. 

Manna Model. BTW model is a deterministic one. Manna pro- 
posed the stochastic sand-pile model jHl by introducing randomness in 
the dynamics of sandpile growth. Here, the critical height is 2. There- 
fore at each toppling the rejected two grains choose their host among the 
four available neighbors randomly with equal probability. After constant 
adding of sand grains, the system ultimately settles at a critical state. 

We consider now the Manna model on a square lattice of size L x L, 
where the sites can be either empty or occupied with unit height i.e., 
the height variables can have binary states hij = 1 or hi j = 0. A site 
is chosen randomly and one height is added at that site. If the site is 
initially empty, it gets occupied: 

hj hj + 1> (1-2) 

If the chosen site is previously occupied then a toppling or 'hard core 
interaction' rejects both the heights from that site: 

hij hij - 2, (1.3) 

and each of these two rejected heights stochastically chooses its host 
among the 4 neighbors of the toppled site. The toppling can happen in 
chains if any chosen neighbor was previously occupied and thus cascades 
are created. After the system attains stable state (dynamics stopped), 
a new site is chosen randomly and unit height is added to it. Thus 
the system evolves in discrete time steps. Here again the boundary is 
assumed to be completely absorbing so that heights can leave the system 
due to the toppling at the boundary. 

With a slow rate of addition of heights (sand grains) at random sites, 
initially the average height of the system grows with time and soon the 
system approaches toward a critical average height he (see Fig. 1(b)). 
Here also the critical average height he has a finite size dependence and 
a similar finite size scaling fit hc{L) = hc{oo) + CL~^^'^ gives v ~ 1.0 
and he = hc{oo) ~ 0.716 (see inset of Fig. 1(b)). This is close to an 
earlier estimate he — 0.71695 [Jj, made in a somewhat different version 
of the model. The avalanche size distribution has got power laws similar 
to the BTW model, at this self-organised critical state at hav = he- 
However the exponents seem to be different ^5^ for this stochastic model, 
compared to those of BTW model. 



Precursors of catastrophic failures 



5 



Precursors of the SOC point: 

In the BTW model. At an average height hav (< hc{L)), when 
ah sites of the system have become stable (dynamics have stopped), a 
fixed number of height units /ip = 4 (pulse of sand grains) is added 
at any central point of the system. Just after this addition, the local 
dynamics starts and it takes a finite time or iterations (r) to return back 
to the stable state (/ijj < 4 for all (i, j)) after several toppling events 
A; and the disturbance spreads over a length ^, the correlation length 
of the system jBl E] . All these three parameters are seen to diverge 
as hav approaches the critical height from below following power laws 
r ~ (/ic - hav)~^, where 7 ^ 1.2; A ~ (/i^ - hav)~\ where 5 ^ 2.0; i ~ 
{he - hav)'", where v ^ 1.0 (see Fig. 2). 
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Figure 2. The variations of the precursors with hav (< hc(L)) in the BTW model for different 
system sizes: L = 100 (plus) L = 200 (cross) and L = 300 (open circle), (a) For relaxation 
time r; in the inset t~^-^ is plotted against hav (b) For the total number of topplings A; 
inset shows A"'^'^ versus hav plot, (c) For the correlation length ^; in the inset, ^"^ "^ is 
plotted against hav 
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In the Manna Model. The precursor parameters in Manna model 
have been measured using similar method as in BTW model. Here the 
pulse of height is hp = 2 and is added to any arbitrary central site. We 
get exactly similar power law behaviors for all the parameters as the 
critical point is approached from below: r ~ {hc — hav)~^-, where 7 = 1.2; 
A ~ (/ic - Kv)~\ where 5 = 2.0; ^ ^ {he - hav)'" , where u = 1.0 (see 
Fig. 3). 
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Figure 3. The variations of the precursors with hav (< hc{L)) in the Manna model for 
different system sizes; L = 100 (plus) L = 200 (cross) and L = 300 (open circle), (a) For 
relaxation time t; in the inset t~^'^ is plotted against hav (b) For the total number of 
topplings A; inset shows A""-^ versus hav plot, (c) For the correlation length ^; in the inset, 
^—1.0 jg plotted against hav 

The Monte Carlo studies showed that these response parameters like 
the relaxation time r, size of the damage A and its radial size ^, all 
tend to diverge, following the above mentioned robust power laws, in 
both the models. Precise knowledge of these power laws can therefore 
help estimating the critical or catastrophic point (he) by extrapolating 
the inverse power of these quantities for informations at hav (< he) (see 
insets of Figs. 2 and 3). Similar critical slowing down phenomena should 
be of interest in the context of magnetohydrodynamics of plasma [Hj. 
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3. 



Precursors of fracture-failure in fiber bundles 



The fiber bundle model was introduced by Peirce :2 in the context of 
testing the strength of cotton yarns. Fiber bundles are of two classes 
with respect to the time dependence of fiber strength: The 'static' bun- 
dles contain fibers whose strengths are independent of time, where as 
the 'dynamic' bundles are assumed to have time dependent elements 
to capture the creep rupture and fatigue behaviors. According to the 
load sharing rule, fiber bundles are being classified into two groups: 
Global load-sharing (democratic) bundles and local load-sharing bun- 
dles. In democratic bundles intact fibers bear the applied load equally 
and in local load-sharing bundles the terminal load of the failed fiber 
is given equally to all the intact neighbors. With steadily increasing 
load, the fiber bundles approach their respective failure point obeying 
a dynamics determined by the load sharing rule. The phase transition 
[01 and dynamic critical behavior of the fracture process in such demo- 
cratic bundles has been established through recursive formulation [21 [HI 
of the failure dynamics. The exact solutions of the recursion relations 
in the global load-sharing case suggest universal values of the exponents 
involved. 

We discuss here the failure of a static fiber bundles under global load 
sharing (democratic bundles) with steadily increasing load on the bun- 
dles. We show analytically the variation of associated precursor param- 
eters with the applied stress which help to estimate the failure point 
accurately. 

The model assumes equal load sharing, i.e., the intact fibers share the 
applied load equally. The strength of each of the fibers in the bundle is 
determined by the stress value (cr^h) it can bear, and beyond which it 
fails. The strength of the fibers are taken from a randomly distributed 
normalised density p{crth) within the interval and 1 such that 



The global load sharing assumption neglects 'local' fiuctuations in stress 
(and its redistribution) and renders the model as a mean-field one. 

Breaking dynamics of the democratic bundles. The breaking 
dynamics starts when an initial stress a (load per fiber) is applied on 
the bundle. The fibers having strength less than a fail instantly. Due 
to this rupture, total number of intact fibers decreases and now these 
intact fibers have to bear the applied load on the bundle. Hence effective 
stress on the fibers increases and this compels some more fibers to break. 




(1.4) 
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These two sequential operations, the stress redistribution and further 
breaking of fibers continue till an equihbrium is reached, where either 
the surviving fibers are strong enough to bear the applied load on the 
bundle or all fibers fail. 

This breaking dynamics can be represented by recursion relations in 
discrete time steps. Let Ut be the fraction of fibers in the initial bundle 
that survive after time step t, where time step indicates the number of 
(elemental) stress redistributions. Then the redistributed load per fiber 
after t time step becomes 

= ^; (1.5) 

and after t + 1 time steps the surviving fraction of fiber is 

Ut+i = I - Piat); (1.6) 



where P{(Tt) is the cumulative probability of corresponding density dis- 
tribution p{ath)'- 

Piat) = r p{ath)dath. (1.7) 

JO 

Using now Eqns. (1.5) and (1.6) we can write the recursion relations 
which show how ut and Ut evolve in discrete time: 

^*+^=l-P(a,)'^° = ^ ^^-^^ 

and 

[/t+i = l-P((7M);C/o = l. (1.9) 



The recursion relations (1.5) and (1.6) represent the basic dynamics 
of failure in global load sharing models. At the equilibrium or steady 
state Ut-\-i = Ut = U* and cr^+i = at = a*. This corresponds to a fixed 
point of the recursive dynamics. Eqn. (1.6) can be solved at the fixed 
point for some particular distribution of p{ath) and these solutions near 
U* (or a*) give the detail features of the failure dynamics of the bundle. 



Precursors of global failure: 

For uniform distribution of fiber strength. We choose first 
the uniform density of fiber strength distribution to solve the recursive 
failure dynamics of democratic bundle. Here, the cumulative probability 
becomes 

/■ft /■ft 

Pi(^t) = / p(,crth)dath = / dath = (Tt- (1-10) 
Jo Jo 
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Therefore Ut follows a simple recursion relation (following Eqn. (1.6)) 
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Figure 4- Three simple models considered here assume (a) uniform, (b) linearly increasing 
and (c) linearly decreasing density p{(Tth) of the fiber strength distribution up to a cutofli' 
strength. 

At the equilibrium state {Ut+i = Ut = f/*), the above relation takes 
a quadratic form of U* : 



U* -U* + cj = {). 



:i.i2) 



The solution is 

U*{a) = \±{a,-af''-a, = \. (1.13) 

Here dc is the critical value of the initial applied stress beyond which the 
bundle fails completely. The solution with (+) sign is the stable one, 
whereas the one with (— ) sign gives unstable solution. The quantity 
U*{a) must be real valued as it has a physical meaning: it is the fraction 
of the original bundle that remains intact under a fixed applied stress a 
when the applied stress lies in the range < cr < Uc. Clearly, U*{ac) = 
1/2 (putting (7 = cJc in Eqn. (1.13)). Therefore the stable solution can 
be written as 

U*{a) = U*{a,) + {a, - a)'/^;a, = ^. (1.14) 

For c7 > (Tc we can not get a real-valued fixed point as the dynamics 
never stops until Ut = 0, when the bundle breaks completely. It may be 
noted that the quantity U* {a) — U* (uc) behaves like an order parameter 
that determines a transition from a state of partial failure {cr < a^) to a 
state of total failure {a > ac) [2]: 

O ^ U*{a) - U*{a,) = {a, - af;P= ^. (1.15) 
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To study the dynamics away from criticality (a — > CTc from below), we 
replace the recursion relation (1.11) by a differential equation 



dU _U^ -U + a 

It ~ u ■ 



(1.16) 



Close to the fixed point we write Ut{cr) = U*{a) + e (where e ^ 0). 
This, following Eq. (1.13), gives 0| 



Ut{a)-U*ia)^eM-t/T), 



(1.17) 



where t = ^ 



l(a,-a)-V2 + i 



. Near the critical point we can write 



T (X {ac- cr) 



'"=2- 



(1.18) 



Therefore the relaxation time diverges following a power-law as a ^ ac 
from below |^. 

One can also consider the breakdown susceptibility x, defined as the 
change of U*{a) due to an infinitesimal increment of the applied stress 
a [3 El 



X 



dU*{a) 



da 



-iar — a) 77 = - 
2^ ' 2 



(1.19) 



from Eqn. (1.13). Hence the susceptibility diverges as the applied stress 
a approaches the critical value ac = \. Such a divergence in x had 
already been observed in the numerical studies f^. 



For linearly increasing distribution of fiber strength. 

the cumulative probability becomes 



crt rcrt 

p{ath)dath = 2 / athdath = c^t ■ 
JO 



Therefore Ut follows a recursion relation (following Eq. (1.6)) 

2 



t+i 



1 



Here, 
(1.20) 

(1.21) 



At the fixed point {Ut+i = Ut = U*), the above recursion relation can 
be represented by a cubic equation of U* 



{U*f-{U*f + a^ = 0. 



(1.22) 
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Solving the above equation we get [H] the value of critical stress a"c = 
a/4/27 which is the strength of the bundle for the above fiber strength 
distribution. Here, the order parameter, susceptibility, relaxation time 
all follow the same power laws (Eqns. (1.15), (1-18) and (1.19)) observed 
for uniform strength distribution. 

For linearly decreasing distribution of fiber strength. In this 
case, the cumulative probability becomes 

P{at) = r p{ath)dath = 2 H (1 - ath)dath = 2at - (1.23) 
Jo Jo 

and Ut follows a recursion relation (following Eqn. (1.6)) 

At the fixed point (C/j+i = Ut = U*), the above recursion relation can 
be represented by a cubic equation of U* 

{U*f-{U*f + 2aU*-a'^ = 0. (1.25) 



Solution of the above equation gives [S] the value of critical stress 
(Tc = 4/27 which is the strength of the bundle. Here also, the precursor 
parameters, the order parameter, susceptibility and relaxation time, all 
follow the same power laws (Eqns. (1.15), (1.18) and (1.19)) observed 
for uniform strength distribution. 

Thus the democratic fiber bundles (for different fiber strength distri- 
butions) show phase transition from a state of partial failure to total 
failure with a well defined precursors (order parameter, susceptibility 
and relaxation time) which show similar power law variation on the way 
as the critical point is approached; characterized by the universal values 
of the associated exponents (a, /3and r]). All the above scaling behaviors 
represented by Eqns. (1.15), (1.18) and (1.19) for stress (a) below the 
global failure stress (cc) of the bundle suggests that a prior knowledge 
of the responses like the fraction of failed fibers, relaxation time etc. can 
be extrapolated following the above power laws to estimate the global 
failure point o"c. 



4. Two-fractal overlap model of earthquake and 
the prediction possibility of large events 

The two-fractal overlap model of earthquake [H] is a very recent model- 
ing attempt. Such models are all based on the observed 'plate tectonics' 
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and the fractal nature of the interface between tectonic plates and earth's 
solid crust. The statistics of overlaps between two fractals is not studied 
much yet, though their knowledge is often required in various physi- 
cal contexts. For example, it has been established recently that since 
the fractured surfaces have got well-characterized self-affine properties 
[TU] , the distribution of the elastic energies released during the slips be- 
tween two fractal surfaces (earthquake events) may follow the overlap 
distribution of two self-similar fractal surfaces Chakrabarti and 

Stinchcombe |5] had shown analytically that for regular fractal overlap 
(Cantor sets and carpets) the contact area distribution follows a simple 
power law decay. The two fractal overlap magnitude changes in time as 
one fractal moves over the other. The overlap (magnitude) time series 
can therefore be studied as a model time series of earthquake avalanche 
dynamics [T^ . 

Here, we study the time (t) variation of contact area (overlap) m{t) 
between two well-characterized fractals having the same fractal dimen- 
sion as one fractal moves over the other with constant velocity. We have 
chosen only very simple fractals: regular and random Cantor sets, gas- 
kets and percolating cluster. We analyse the time series data of Cantor 
set overlaps only to find the prediction possibility of large events (oc- 
currence of large overlaps). We show that the time series m{t) obtained 
by moving one fractal uniformly over the other (with periodic boundary 
condition) has some features which can be utilized to predict the 'large 
events'. This is shown utilizing the discrete or a 'quantum' nature of the 
integrated (cumulative) overlap over time. 

Two-fractal overlaps: 

We now discuss the two- fractal overlap sizes. Let us take two identical 
Cantor sets at finite generation and slide one over the other, assuming 
periodic boundary condition. Two such cases for a regular Cantor sets 
as well as a gasket are shown in Fig. 5. Similar cases for random fractals 
are shown in Fig. 6. 

(a) (b) 

rz IT 

n=3 n=3 
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Figure 5. (a) A regular Cantor set of dimension In 2/ In 3; only three finite generations are 
shown, (b) The overlap of two identical (regular) Cantor sets, at n = 3, when one slips 
over other; the overlap sets are indicated within the vertical lines, where periodic boundary 
condition has been used, (c) A regular gasket of dimension In 3/ In 2 at the 7th generation, 
(d) The overlap of two identical regular gaskets at same generations (n = 7) is shown as one 
is translated over the other; periodic boundary condition has been used for the translated 
gasket. 
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Figure 6. (a) A random Cantor set of dimension In 2/ In 3; only three finite generations are 
shown, (b) Overlap of two random Cantor sets (at n = 3; having the same fractal dimension) 
in two different realisations. The overlap sets are indicated within the vertical bars, (c) A 
random realisation of a gasket of dimension In 3/ In 2 at 7th generation, (d) The overlap of 
two random gaskets of same dimension and of same generation but generated in different 
realisations. 



Next, we study the overlap distribution of two well-characterised random 
fractals; namely the percolating fractals ^3]. It seems, although many 
detailed features of the clusters will change with the changes in (parent) 
fractals, the subtle features of the overlap distribution function remains 
unchanged. 
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Figure 7. The overlap between two percolating clusters; (a) and (b) are two typical realisa- 
tions of the same percolating fractal on square lattice (rf/ ~ 1.89) and (c) their overlap set. 
Note, the overlap set need not be a connected one. 

Our earlier study [H] on the overlap statistics for regular and random 
Cantor sets, gaskets and percolating clusters indicated a universal scaling 
behavior of the overlap or contact area (m) distributions P{m) for all 
types of fractal set overlaps mentioned: P'{m') = L"P{m, L);m' = 
mL~"', where L denotes the finite size of the fractal and the exponent 
a = 2{df — d); df being the mass dimension of the fractal and d is the 
embedding dimension. Also the overlap distribution P{m), and hence 
the scaled distribution P'{m'), are seen to decay with m or m' following 
a power law (with exponent value equal to the embedding dimension 
of the fractals) for both regular and random Cantor sets and gaskets: 
P[m) ~ m~^;P = d. However, for the percolating clusters the 
overlap size distribution takes a Gaussian form. 

We consider now the time series obtained by counting the overlaps 
m{t) as a function of time as one Cantor set moves over the other (peri- 
odic boundary condition is assumed) with uniform velocity. The analysis 
presented here follows our previous study 
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Overlap time series data. The time series are shown in Fig. 8, 
for finite generations of Cantor sets of dimensions In 2/ In 3 and In 4/ In 5 
respectively. 
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Figure 8. The time (t) scries data ofoverlap size (m) for regular Cantor sots: (a) 
of dimension In2/ln3, at 8th generation; (b) of dimension In4/ln5, at 6th genera- 
tion. The obvious periodic repeat of the time series comes from the assumed periodic 
boundary condition of one of the sets (over which the other one shdes) . 



Cumulative overlap sizes. Now we calculate the cumulative 

overlap size Q{t) = J^mdt and plot that against time in Fig. 9. Note, 
that 'on average' Q(t) is seen to grow linearly with time t for regular as 
well as random Cantor sets. This gives a clue that instead of looking 
at the individual overlaps m{t) series one may look for the cumulative 
quantity. We observe Q{t) ~ ct, where c is dependent on the fractal. 
This result is even more prominent in the case of Cantor sets with df = 
In 4/ In 5. 
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Time (t) Time (t) 

Figure 9. The cumulative overlap Q versus time; for pure cantor sets: (a) of dimension 

In 2/ In 3 (at 8th generation) and (b) of dimension In 4/ In 5 (at 6th generation). The 
dotted line corresponds to those for two identical but random Cantor sets. In (b) the 
two lines fall on each other. 
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Cumulative overlap quantization. We first identify the 'large 

events' occurring at time tj in tlie m{t) series, where m{ti) > M, a 
pre-assigned magnitude. Wc then look for the cumulative overlap size 
Qi{t) = J^._^ mdt, t < ti, where the successive large events occur at times 
ti-i and U. The behavior of Qi with time is shown in Fig. 10 for regular 
cantor sets with df = In 2/ In 3 at generation n = 8. Similar results 
are also given for Cantor sets with df = In 4/ In 5 at generation n = 6 
in Fig. 11. Qi{t) is seen to grow almost linearly in time up to Qi{ti) 
after which it drops down to zero. It appears that there are discrete 
(quantum) values of Qi{ti). One can therefore anticipate large events 
when the cumulative seismic overlap Qi{t), since the last major event 
at ti-i, grows linearly with time t to the discrete levels IQo, specific to 
the series of events. No major event occurs during the growth period of 
Qi{t) for the intermediate values. 




1000 2m 3000 4000 5000 6000 70OO 1000 2TO 3000 4000 5000 6000 7000 
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Figure 10. The cumulative overlap size variation with time (for regular Cantor sets of 
dimension In 2/ In 3, at Sth generation), where the cumulative overlap has been reset 
to value after every big event (of overlap size > M where M = 128 (a) and 32 (b) 
respectively). 
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Figure 11. The cumulative overlap size variation with time (for regular Cantor sets of 

dimension In 4/ In 5, at 6th generation), where the cumulative overlap has been reset 
to value after every big event (of overlap size > M where M = 2400 (a) and 2048 
(b) respectively). 



5. Summary and Discussions 

In all the dynamical models of sandpile avalanches (discussed in section 
2), we find that the growing correlations in the dynamics of constituent 
elements manifest themselves as various precursors. The number of top- 
plings A, relaxation time r and the correlation length i^, in both BTW 
and Manna model, grow and diverge following power laws as the sys- 
tems approach their respective critical points he from below (see Figs. 
2 and 3): A ~ (/ic - hav)~^ , t ^ {he - hav)~"' and ^ ~ (/ic - hav)'" ■ 
For two dimensional systems, we find numerically here S ~ 2.0, 7 ~ 1.2 
and u ~ 1.0 for both BTW and Manna model. Basically, we studied 
the behavior for hav < he, the precursor behavior, where ^ is necessarily 
finite. As we add here the tiny pulse at some central site of a relatively 
large system, the boundary effect can not be really felt because of the 
smallness of ^ compared to L for most values of hav This explains the 
lack of finite size effect in our precursor studies. 

For the global failure or fracture in fiber bundle models (discussed 
in section 3), we find that the breaJcdown susceptibility x (giving the 
increment in the number of broken fibers for an infinitesimal increment 
of load on the bundle) and the corresponding relaxation time r (required 
for the bundle to stabilise, after successive failures of the fibers), both 
diverge as the external load or stress approaches its global failure point 
CTc from below: r ~ (tic -fJai,)"" and x ~ (o"c -CTav)'"^] a = r] = 1/2. 
These results are of course analytically derived for different fiber strength 
distributions and the universality class of the model has been confirmed. 
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If one Cantor set moves uniformly over another, the overlap between 
the two fractals change quasi-randomly with time (see eg., Fig. 8). These 
time series are taken as model seismic activity variations in such two- 
fractal overlap model of earthquake. The overlap size distribution was 
argued 0] and shown to follow power law decay. We showed numeri- 
cally (see section 4) that if one fixes a magnitude M of the overlap sizes 
m, so that overlaps with m > M are called 'events' (or earthquake), then 
the cumulative overlap Qi{t) = ^ mdt, t < ti, (where two successive 
events oim> M occur at times and ti) grows linearly with time up 
to some discrete quanta Qi{ti) = IQq, where Qq is the minimal overlap 
quantum, dependent on M and n (the generation number). Here I is an 
integer (see Figs. 10, 11) ^3]. Although our results here are for regular 
fractals of finite generation n, the observed discretisation of the overlap 
cumulant Qi with the time limit set by n, is a robust feature and can be 
seen for larger time series for larger generation number n. This part of 
the model study therefore indicates that one can note the growth of the 
cumulant seismic response Qi{t), rather than the seismic event strength 
m{t), and anticipate some big events as the response reaches the discrete 
levels IQo, specific to the series of events. 

Some prior knowledge of the precursors and their precise behavior 
(power laws etc.), as discussed here for three different kinds of failure 
models of cooperatively interacting dynamical systems which are prone 
to major failures, are therefore plausible. These precursor responses 
should help estimating the location of the global failure or critical point 
using proper extrapolation of the above quantities, which are available 
long before the failure occurs. The usefulness of such precursors can 
hardly be overemphasized. 
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